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Abstract 

Mostly for the purpose of applications for the energy and the environment and 
for the design of sources of neutrons or exotic nuclides, intense research has been 
dedicated to spallation, induced by protons or light projectiles at incident energies 
of around 1 GeV. In this energy range, while multifragmentation has still a minor 
share in the total reaction cross section, it was observed to have, together with fis- 
sion, a prominent role in the production and the kinematics of intermediate-mass 
fragments, so as to condition the whole production of light and heavy nuclides. The 
experimental observables we dispose of attribute rather elusive properties to the 
intermediate-mass fragments and do not allow to classify them within one exclu- 
sive picture which is either multifragmentation or fission. Indeed, these two decay 
mechanisms, driven by different kinds of instabilities, exhibit behaviours which are 
closely comparable. High-resolution measurements of the reaction kinematics trace 
the way for probing finer features of the reaction kinematics. 
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1 Introduction 



Sixty years passed from Serber's early description [1] of high-energy nuclear 
reactions induced by nucleons and light nuclei. In their main outline, such re- 
actions are describes as the exciting of an atomic nucleus, followed by a decay 
process in several nuclides, clusters, protons and neutrons [2] . Several different 
nuclear systems have been explored with beams of various incident energies, 
at different facilities and with different experimental techniques. Numerous 
experimental results and various observables on the production of residues, on 



the kinematics of the emission of ejectiles and fragments, and on the correla- 
tions inspired and constrained the physical models. 

At incident energies of few hundred MeV per nucleon, an excited and fully 
equilibrated complex, named compound nucleus, is formed and successively 
de-excites by mainly fission-evaporation decays. When the excitation energy 
of the hot nucleus exceeds the threshold for emission of particles or clusters 
(including fission), the system has the possibility to decay by any of the open 
channels. If the excited system is not too hot, the favoured process is the re- 
ordering of its configurations: a great number of arrangements are available 
where all nucleons remain in states below the continuum, occupying excited 
single-particle levels around the Fermi surface. Rather seldom, compared with 
this thermal chaotic motion of the system, one nucleon acquires enough energy 
to pass above the continuum and may eventually leave the nucleus. This pic- 
ture was extended to include the production of intermediate-mass fragments 
by cluster decay and oscillations in fission direction as well. In this respect, 
there would be a gradual transition from very asymmetric to symmetric con- 
figurations in the binary split of the decaying compound nucleus, so that evap- 
oration of nucleons and light nuclei on the one hand and symmetric fission on 
the other hand are just the opposite extremes of the manifestation of the same 
process. This generalisation, introduced by Moretto [3,4], allows to name fis- 
sion in a generalised sense all (binary) decays of a compound nucleus. Since 
this decay is a rare process, one evaporation event, or fission event, proceeds 
after the other, sequentially. 

At incident energies of some GeV per nucleon, a very highly excited compos- 
ite nuclear system is formed; it undergoes a violent de-excitation processes, 
named multifragmentation, which manifests by the production of several mas- 
sive fragments filling broad kinetic-energy spectra. General reviews on this 
process can be found in refs. [5,6]; a more specific review treating multifrag- 
mentation induced with high-energy proton beams can be found in ref. [7]. The 
amount of experimental evidence suggests to describe the de-excitation as a 
simultaneous disassembling of the hot nucleus in several constituents. From 
the experimental observation of angle correlations between fragments, it was 
found that the emission of two fragments with small relative angles is highly 
improbable. Such a dependence on relative emission angles is on the contrary 
absent for the ejectiles of a sequential evaporation process [8]. The disintegra- 
tion is so rapid to be considered "simultaneous" , in the sense that it evolves 
in so short a time interval (10 -22 -10~ 21 s) that fragments can still exchange 
interactions while they are accelerated in their mutual Coulomb field. On the 
other hand, the process is expected to be still sufficiently "slow" to exceed the 
relaxation time of the strong interactions and, for this reason, the system is 
assumed to be thermalised before it disintegrates. Within a thermodynamical 
picture, multifragmentation is a phenomenon related to the equation of state 
of nuclear matter [9,10]; with a certain resemblance with Van-der-Waals fluid, 
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hot bulk matter would enter the liquid-gas coexistence region of the phase di- 
agram and separate in the corresponding coexisting liquid and gas phases (see 
refs. [12,11] for review), driven by a rapid amplification of spinodal instabili- 
ties [13]. According to the same picture, when the composite system reaches 
the conditions for disintegrating, it should be diluted to some fractions of the 
nuclear saturation density. These low densities are reached as a result of a dy- 
namical process of expansion, which is expected to explain the high velocity 
of the fragments observed experimentally. There are however also alternative 
mechanisms proposed to explain the experimental features of the multifrag- 
mentation process which are close to the fully equilibrated compound-nucleus 
decay [14,15]. 

In recent years, mostly for the purpose of energetic and environmental appli- 
cations, as well as for the production of beams of neutrons or exotic nuclides 
in accelerator facilities, increasing interest is devoted to spallation reactions 
induced by protons or deuterons at incident energies close to 1 GeV. At this 
energy, the reaction is situated somehow in between the two scenarios de- 
scribed above and is particularly interesting for studying the transition from 
fission to multifragmentation. The physics underlying the spallation process 
in this energy range is however rather uncontrolled. From the experimental 
side, the approaches used up to now have always provided a partial survey 
of all the observables which are necessary to formulate conclusive answers. 
Inclusive approaches are best suited for the identification of the nuclides, for 
the measurement of the corresponding production cross sections and kinetic- 
energy distributions; they however neither provide particle correlations, like 
the multiplicity of the fragments and particles formed in the reaction, nor an- 
gle (or velocity) correlations, necessary to probe the reaction kinematics for 
each single event. Correlation observables are the specificity of exclusive ap- 
proaches, based on the employment of multidetectors which, at the moment, 
are still not able to provide the full isotopic identification up to heavy ele- 
ments and high-resolution velocity measurements. On the other hand, from 
the modelling side, several different strategies based either on a generalised 
fission-evaporation scenario [3,16] or on a statistical description of multifrag- 
mentation [17,18,19,20,21,22,23,24] have been able to describe the same ex- 
perimental data with comparable accuracy. 



2 Residue production in spallation induced by 1 GeV protons 

The measurement of the production of spallation residues at relativistic en- 
ergy has been the purpose of several years of experiments performed at the 
FRagment Separator [25,26] (GSI, Darmstadt). The reactions were measured 
in inverse kinematics, by directing heavy-ion beams at various energies on a 
target of liquid hydrogen or deuterium. Fig. 1 presents a survey on the nuclide 
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production for the spallation reactions induced by protons at 1 GeV which 
were measured at the FRagment Separator. 

For the reaction 238 U(i J 4GeV) +P [27] the nuclide distribution reflects the in- 
terplay of more decay patterns. Quite independently on the entrance chan- 
nel (and in particular the neutron enrichment of the compound nucleus), 
evaporation residues [28] populate the neutron-poor side of the nuclide chart 
around the residue corridor [29], which corresponds to the situation where [30] 




Fig. 1. Experimental survey on the distribution of spallation residues pre- 
sented in nuclide charts for the systems UriAGeV) + P [27,28,31,32,33], 
208 Pb ( iAGeV) +P [34,35], 136 Xe( 1AGeV) +p [40], and 56 Fe (1AGeV) +P [37,36]. Filled 
circles indicate the projectiles (in inverse kinematics). 
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dN/dZ = (Tn/Tz), where Tz and are the proton- and neutron-emission 
width, respectively. This is for instance the reason why neutron-rich nuclei 
are favoured materials for neutron sources. In the centre of the nuclide distri- 
bution an imposing high-energy fission hump emerges between two shoulders 
generated by low-energy fission channels [31], which are exotic and asymmet- 
ric. Other asymmetric fission channels are determined by high-energy fission, 
as an extension of the central hump; they overlap with the evaporation pro- 
duction on one side [32] and extend to the intermediate-mass fragments [33] 
on the opposite side. 

While uranium is an attractive material for the production of exotic nuclei 
by low-energy fission, lead is preferred as a neutron source because, as the 
measurement of 208 Pb(iAGeV) +P reveals [34,35], fission is much reduced with 
respect to evaporation channels. 

Iron was measured because it is a structural material in nuclear installations. 
From a phenomenology point of view, the results on the reaction 56 Fe(i AGeV) + 
p [37,36] brought the attention to light systems, below the Businaro-Gallone 
point [38,39]. At variance with heavy systems, the fission potential becomes 
convex and fission channels are asymmetric. Moreover, the excitation energy is 
larger in lighter systems when reactions are compared at the same incident en- 
ergy. As a result of reaction models [37], for small impact parameters the exci- 
tation energy of the composite system generated in the reaction 56 Fe(i AGeV)+P 
is found to even exceed 3 MeV per nucleon; around this value multifragmenta- 
tion is expected to set in. For the production of intermediate-mass fragments, 
this process sums up to the asymmetric fission channels. 

The measurement of 136 Xe(i agcV) +P was performed on purpose to analyse the 
intermediate-mass-fragment production and its influence on the overall decay 
mechanism. In this respect, 136 Xe offers an optimum point of observation: sym- 
metric fission, which in heavier systems hides the intermediate-mass-fragment 
production, is suppressed because the system is slightly below the Businaro- 
Gallone point; then, below the Businaro-Gallone point, 136 Xe is the stable 
nuclide with the largest neutron excess N — Z, and it approaches the neutron 
enrichment N/Z of 208 Pb, so that a comparison can be made between the two 
systems; and finally, the system 136 Xe(i agcV)+P is excited right enough to still 
allow for some oscillations towards multifragmentation. The results on the nu- 
clide production [40] illustrate that, besides a close similarity with the system 
208 Pb( 1 ^ GcV ) +p for the evaporation features, beyond a mass loss of around 
A A = 70, the ridge of the residue production abandons the neutron-poor side 
of the nuclide chart around Z = 40 and migrates progressively towards the 
neutron-rich side for lighter residues. The lightest residues even populate the 
neutron rich side of the nuclide chart with respect to the valley of stability. 

It may be remarked that both asymmetric fission channels and multifragmen- 
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tation of a neutron-rich composite system produce residues with high neutron 
excess in the average; the successive evaporation of residues, which is expected 
due to the high mean excitation energy of the system, contributes to further 
dissipate the traces of the initial stage of the decay process [41]. Qualitatively, 
both the two mechanisms also result in a similar U-shape of the mass distri- 
bution of residues, with decreasing depth of the hollow for increasing mean 
excitation energy of the system. For this reason, the residue production alone 
is a rather elusive observable and it could be reproduced with comparable 
quality both on the basis of multifragmentation processes and asymmetric 
fission. 



3 Kinematics 

More robust observables for the reaction mechanism are those related to the 
Coulomb field experienced by the fragments. In exclusive experiments, these 
are angle correlations among fragments, which reflect the role of the Coulomb 
field in each event; such experimental strategy, based on multidetectors and 
largely used for ion-ion collisions or for multifragmentation induced by high- 
energy protons, was unfortunately not employed in the study of spallation 
in the incident-energy range of around 1 GeV. In inclusive experiments, the 
observables probing the Coulomb field are the mean quantities related to the 
momentum distribution of fragments. Evidently, inclusive observables mixes 
up the contributions of all processes responsible for the formation of one given 
fragment; therefore, the identification of the reaction process can not be im- 
posed from a selection of separate observables (like the fragment multiplicity 
or the transfered energy of light charged particle in exclusive experiments [42]), 
but it should be extracted by the momentum distribution itself. Up to a cer- 
tain extent, this is possible with magnetic spectrometers operated in inverse 
kinematics, as in the case of the experiments discussed in section 2, where the 
high-resolution momentum distributions could be deduced from the measure- 
ment of the magnetic rigidity Bp, as precise as 5- 10~ 4 (FWHM) for individual 
reaction products. 

A dedicated analysis procedure (refs. [37,40] for detailed discussion) was ap- 
plied to reconstruct invariant cross sections from the inclusive measurement of 
the momentum distribution at the FRagment Separator. An overview of the 
result is shown in fig. 2 for some light resides of the reaction 136 Xe(iAGeV) +P- 
More details of this observable are shown for the nuclide 20 F: in the insert 
(a), a planar cut along the beam axis of the full velocity distribution in the 
projectile frame is shown in the v± x v\\ space (v± and v\\ are respectively 
the perpendicular and parallel velocity components of the fragment 20 F in the 
beam frame), as reconstructed from the measured momentum distribution; by 
selecting the velocities aligned along the beam axis, this representation can 
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be reduced to the distribution of invariant cross section o\ as a function of v\\ 
in the projectile frame, shown in the insert (b). The distribution of invariant 
cross section <j\ of all the intermediate mas fragments formed in the reaction 
136 Xe(iAGeV) + V results from the overlap of two shapes: a component with a 
convex centre and a component with a concave centre. 



A "convex" mode describes two completely different situations. In a first case, 




2 
Uboost [cm/ns] 



Fig. 2. Distribution of invariant cross section for some intermediate mass fragments 
produced in the reaction 136 Xe(! AGeV) + P [40]- The corresponding nuclides are 
indicated on a portion of the fragment distribution, on the nuclide chart. Two 
kinematical modes, one convex in the centre, the other concave in the centre, are 
indicated. For the nuclide 20 F more details are shown: (a) The planar cut along 
the beam axis of the reconstructed full distribution da/dv h in the beam frame as a 
cluster plot. The concave mode is indicated by the corresponding circular ridge of 
radius w pea k- (b) The invariant-cross-section distribution, as for the other nuclides. 
(C) Reconstructed cross-section distribution for the concave mode as a function of 
the boost velocity in the source frame v boost- The plot allows to calculate the mean 
boost (vboost)j which differs from u pe ak- 



7 



it describes the velocity distribution of evaporation residues after sequential 
emission of nucleons and clusters. At these incident energies, the most prob- 
able excitation energy is just sufficient for the emission of few nucleons and 
the probability of higher excitations, connected with longer evaporation paths, 
can only decrease progressively; this behaviour ensures that the distribution 
of (heaviest) evaporation residues drops in cross section monotonically when 
moving away from the projectile and its contribution to the production of 
intermediate mas fragments is invisible for all the systems described in fig. 1. 
For these systems, a convex distribution of invariant cross section associated 
to intermediate mass fragments rather indicates a multifragmentation pro- 
cess; in particular, the convex mode probes the Coulomb field produced by 
the disassembling of the system in more fragments having a comparable size. 
The explanation for the asymmetry of this mode, characterised by a long tail 
in backward direction and a maximum shifted forward, may be searched in 
possible dynamical effets. It may be interesting to remark that the convex 
mode, associated to multifragmentation, can not be observed in direct kine- 
matics because a threshold in kinetic energy hides the centre of the distribu- 
tion of invariant cross section. The "concave mode" reflects a strong Coulomb 
boost as experienced by fission fragments. When this pattern is associated to 
intermediate-mass fragments, it may indicate an asymmetric fission process as 
well as the possible extension of multifragmentation events to very asymmet- 
ric partitioning configuration; in this latter case, the kinematics is determined 
by the repulsion exerced by one very heavy fragment and it should not be 
disturbed greatly by the presence of more than one light fragment. In con- 
clusion, a concave component in the distribution of invariant cross section 
of intermediate-mass fragments indicates that the system certainly broke up 
asymmetrically, but the attribution of the process to fission or multifragmen- 
tation stays uncertain on the basis of this inclusive observable. 

In fig. 2, the insert (c) represents for the nuclide 20 F the convex component 
alone in polar coordinates in the frame of the corresponding average emitting 
source. The mean value of this spectrum (wboost) would be consistent with a 
fission barrier if the process were exclusively reduced to fission. Such test is 
shown in fig. 3 for some light elements produced in the system 136 Xe(i AGeV)+P- 
The figure compares the measured mean boost (fboost) with the highest ex- 
pected value, determined by the split of the heaviest possible system, 136 Xe. 
The expected boost is the highest when the fission fragments are considered 
not deformed and subjected to the fusion potential, as described by the em- 
pirical formula of Bass [43,?]. The lowest value is calculated when the fission 
fragments are considered deformed and joined through a neck according to the 
liquid-drop picture (scission-point model of Wilkins et al. [45,?]). In between, 
is situated the value calculated from a systematics of total kinetic energy for 
light fissioning nuclei [47], further modified in ref. [37] to describe asymmet- 
ric splits. The lighter are the elements, the more all the three expectations 
underpredict the data. This test suggests that the concave mode should be 
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alimented, in addition to fission, also by multifragmentation channels which 
could be related to an expansion process and could enhance the mean boost. 

Similar experimental observations were collected for the system 56 Fe(iAGeV) + 
p [37]. In the system 238 U(iAGeV) + P, the intermediate-mass fragments were 
fully described by the concave mode; also in this case, the corresponding mean 
boost was characterised by very high values [331 so that it is tempting to make 
a comparison with the lighter systems 56 Fe ( !i eV) + p and -Xe (1 / GeV) +p. 



4 Concluding remarks 

In conclusion, the coupling of kinematical and production observables is a 
rather robust approach to study the reaction mechanism in such a complex 
situation as spallation reactions at incident energies of around 1 GeV. This 
approach indicated that multifragmentation is a possible decay pattern and 
determines the kinematics by producing a convex component in the distribu- 
tion of invariant cross section. More difficult is to define its possible extension 
to the asymmetric breakups of the system, determined by a concave compo- 
nent in the distribution of invariant cross section. Is the observation of a very 
high mean boost in the asymmetric splits a signature of dynamical effects, like 
expansion, which goes mostly in the direction of multifragmentation? 

The left panel of fig. 4 shows the mass distribution of the residues of the system 
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Fig. 3. Evolution of the reconstructed mean boost, deduced for the convex kinetic 
mode as a function of the neutron number for some light elements in the reaction 
136 Xe( 1 ^ GeV ) +p [40] . The measurement is compared with three expectations for the 
Coulomb boost in the split of 136 Xe, according to the total-kinetic-energy system- 
atics of Tavares and Terranova [47], the scission-point model of Wilkins et al. [45,?] 
and the nucleus-nucleus fusion model of Bass [43,?]. 



9 



100 



10 



<J (A ) [mb] 



0.1 




Xe(iAGeV)+P f! 




J i I i L 



ft* 

/ /cf 

f '< 

J ! 



100 



10 




ct(AA) [mb] 



\ \ 208 

w 



V 



J I I I I 



0.1 



Pb(lAGeV)+P | 



°Xe ( iAGeV)+P 

J i I i_Ll U i I i I i I 



20 40 60 80 100 120 140 20 40 60 80 100 120 140 
A AA 

Concave mode 
Convex mode 



Fig. 4. Left panel. Projections of the residue-production cross sections along the 
mass number for the system 136 Xe(i^GeV) +P [40]. For intermediate mass frag- 
ments the contributions from the two kinematical modes, concave and convex, are 
indicated. A Weisskopf calculation for the expected contribution of the evaporation 
of protons, neutrons and alpha particles is shown. Right panel. Mass-loss distribu- 
tion of production cross sections. The experimental results obtained for the systems 
136 Xe(! AGeV) +P [40] and 208 Pb( lv 4GeV) +P [34] are compared. The same Weisskopf 
calculation for the system 136 Xe(! AGeV) +P shown in the left panel, is repeated for 
comparison. 



136 Xe(iAGeV) +P- It a l so shows for the concave and convex kinematical modes 
the corresponding contribution to the production cross section. The limitations 
of the experimental approach did not allow to complete the distribution for 
these two components over the whole mass distribution. However, we can 
appreciate a rise in cross section of the convex mode up to joining the value of 
the concave mode. The convex mode is expected to extend further and to be 
the dominant mode in the region of the hollow of the mass distribution. The 
concave mode, which gives the highest contribution in the production of the 
light nuclides, drops steeply in cross section in proximity of the hollow. Since 
this mode supposes the existence of a heavy partner in a split of mostly binary 
kind, the corresponding distribution of cross section should also aliment the 
region of heavy masses. The side of the mass spectrum in the region of heavy 
masses is often taken as a reference for testing the modelling of the initial 
stage of the reaction process (cascade), due to its sensitivity to the excitation 
energy introduced in the system during the collision. However, already at 
incident energies of around 1 GeV, the contribution of the convex mode to 
the production of heavy masses should greatly modify the slope of the mass 
distribution for the heavy-fragment side of the hollow, with respect to the 
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function characterising a more simple mechanism where only protons, neutrons 
and alpha particles are emitted. Such function, which can be estimated by a 
Weisskopf calculation, drops to imperceptible cross sections for a mass loss 
AA of less than sixty units and deviates form the experimental distribution 
already after A A ps 40. The right panel of fig. 4 shows that the cross section 
evolves with the mass loss in a very similar way for the system 136 Xe(i AGeV)+P 
and 208 Pb(iAGcV) + £> in the region of heavy masses. 

This feature reinforces the idea that the concave mode, associated to a fast 
fission-like process, and observed in the velocity spectra for systems of different 
size, from 56 Fe(iAGeV) + P to 238 UnAGeV) + £>, is a rather general picture in 
spallation reactions induced by 1 GeV protons. On the other hand, the concave 
mode, more evidently related to multifragmentation, is a characteristic of more 
excited system. 

The purpose of this report was to outline the main phenomenological features 
of spallation in the incident-energy range of 1 GeV per nucleon, without dis- 
cussing any consequence on applications. It however derives naturally that the 
passage from the phenomenology to the modelling is necessary for application 
purposes, in order to better describe the overall production of nuclides and 
the kinematics of light fragments in spallation reactions which are and will be 
largely employed in energetic and environmental applications. 
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